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Abstract 



We study the effective second-harmonic generation (SHG) susceptibihty 
in graded metahic films by invoking the local field effects exactly, and further 
numerically demonstrate that the graded metallic films can serve as a novel 
optical material for producing a broad structure in both the linear and SHG 
response and an enhancement in the SHG signal. 
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There are much recent interest and practical need for nonhnear optical materials that 
process large nonlinear susceptibility or optimal figure of merit. A large enhancement in non- 
linear responses has been found for a sub- wavelength multilayer (i.e., thin film) of titanium 
dioxide and conjugated polymer [1]. For nonhnear effects other than the Kerr effect, Hui 
et al. [2] derived general expressions for the effective susceptibility for the second-harmonic 
generation (SHG) in a binary composite of random dielectrics. They have also studied the 
thickness dependence of effective SHG susceptibility in films of random dielectrics and in 
composites with coated small particles [3,4]. Graded materials with various functionalities 
appear in nature and in fabricated materials. Graded thin films have many applications as 
the gradation profile provides an additional control on the physical properties. Graded thin 
films often possess different optical properties [5], when compared to bulk materials. It is 
known that graded materials have quite different physical properties from the homogeneous 
materials [6]. Also, it has been observed that the compositionally graded barium strontium 
titanate thin films have better electric properties than a single-layer barium strontium ti- 
tanate film with the same compositions [7]. How to achieve enhanced SHG is up to now 
a challenging agenda [8,9]. In this paper, we shall investigate SHG in a graded metallic 
film with an intrinsic SHG response and a graded linear response in the metallic dielectric 
function. 

Let us consider a graded metallic film with thickness L, with the gradation profile in the 
direction perpendicular to the film. If we only include quadratic nonlinearities, the local 
constitutive relation between the displacement field 'D{z) and the electric field E(2;) in the 
static case would be Di[z) — J2j ^ij{z)Ej{z) + J2jkXijk{z)Ej{z)Ek{z) [3,4] with i — x,y,z, 
where Di{z) and Ei{z) are the ith component of Ii{z) and £(2;), respectively, and Xijk is 
the SHG susceptibility. Here eij{z) denotes the linear dielectric function, which we assume 
for simplicity to be isotropic eij{z) — e{z)Sij. Both e{z) and Xijki^) functions of z and 
describe the gradation profiles. In general, when a monochromatic external field is applied, 
the nonlinearity will generate local potentials and fields at all harmonic frequency. For a 
finite frequency external electric field of the form £'0 = Eo{u!)e~^'^^ + c.c, the effective SHG 
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susceptibility X2uj can be extracted by considering the volume average of the displacement 
field at the frequency 2uj in the inhomogeneous medium [2-4]. Next, we adopt a graded 
dielectric profile that follows the Drude form 

f ^ 1 ^p(^) /IN 

eiz,u!) = l -. — ^ . , , , , (1) 

where ^ < z < L. The general form in Eq. (1) allows for the possibility of a grada- 
tion profile in the plasma frequency [e.g., Eq. (4)] and the relaxation rate [e.g., Eq. (5)]. 
For a z-dependent profile, we can make use of the equivalent capacitance for capacitors 
in series to evaluate the effective perpendicular linear response for a given frequency, 
l/e(a;) = (1/^) ^z\Vl e{z^Lo)\. Using the continuity of the electrical displacement field, 
the local electric fields E{z^ iS) satisfies 

e{z,uj)E{z,uj) = l{uj)E^{uj), (2) 

where Eq{cu) is the applied field along z axis. A z-dependent profile for the plasma frequency 
and the relaxation time can be achieved experimentally. One possible way may be to impose 
a temperature profile, as it has been observed that surface enhanced Raman scattering 
is sensitive to temperature [10]. Thus, one may tune the surface plasmon frequency by 
imposing an appropriate temperature gradient [11]. A temperature gradient may also be 
used in materials with a small band gap or with a profile on dopant concentrations. In this 
case, one may impose a charge carrier concentration profile to a certain extent. This effect, 
when coupled with materials with a significant intrinsic nonlinear susceptibility, will give us 
with a way to control the effective nonlinear response. For less conducting materials, one 
may replace the Drude form of dielectric constants by a Lorentz oscillator form. It may 
also be possible to fabricate dirty metal films in which the degree of disorder varies in the 
2;-direction and hence leads to a relaxation-rate gradation profile. 

The calculation of the effective nonlinear optical response then proceed by applying the 
expressions derived in Refs. [2,3]. Next, the effective SHG susceptibility X2uj is given by 
X2u> = {X2uj{z)Eiin{z,2u})Eiin{z,u})^)/[Eo{2u})Eo{u}y], where Eun is the linear local electric 



field in the graded film with the same gradation profile but with a vanishing nonlinear 
response at the frequency concerned. Using Eq. (2) for the linear local fields, the effective 
SHG susceptibility can be rewritten as an integral over the film as 

To illustrate the SHG in graded films, we consider as a model system that the intrinsic 
SHG susceptibility X2uj{z) — Xi to be a real and positive frequency- independent constant and 
does not have a gradation profile. In doing so, we are allowed to focus on the enhancement 
of the SHG response when compared to Xi. To show the effects of gradation, here we take 
as a model plasma-frequency gradation profile 

Up{z)^Upmi-C^-z), (4) 

and a model relaxation-rate gradation profile [12] 

7(z) = 7(00) + Cjz, (5) 

where and are constant parameters tuning the profile which is assumed to be linear in 
z. Here 7(00) denotes the bulk damping coefficient, i.e., for 2; — > 00. Set thickness L = 1 so 
that we could focus on the film with a fixed thickness. Regarding the thickness dependence, 
we refer the reader to the work by Hui et al. [3] 

Figure 1 shows the real and imaginary parts of the effective linear dielectric constant 
[Fig. 1(a) and Fig. 1(b)], and the real and imaginary parts of the effective SHG susceptibility 
[Fig. 1(c) and Fig. 1(d)] as a function of frequency uj/ujp{Q). Also shown is the modulus of 
X2w/Xi) see Fig. 1(e). The dielectric function gradation profile is given in Eqs. (1), (4) and (5) 
with = 0, i.e., only a graded plasmon frequency is included. Throughout the calculations, 
the real part of the (linear) dielectric constant is negative naturally. In this case, a broad 
resonant plasmon band is observed. Note that for C^, — > 0, Lijp{z) / u;p{0) — > 1. As C^, 
increases, ujp{z) takes on values within a broader range across the film, and leads to a broad 
plasmon band. Increasing C^j also causes the plasmon peak to shift to lower frequencies. 
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The reason is that, in analogous to capacitors in series, the effective dielectric constant of the 
film is dominated by the layer with the smallest dielectric constant. For the SHG response, 
the frequency dependence is quite complicated. As C,^ increases, it is noted that structures 
in the SHG response also shifts to the lower frequencies and the range of values of the SHG 
susceptibility increases as well. Figure 2 displays the results of model calculations in which a 
gradation profile of the relaxation rate of the form in Eq. (5) is also included. The effects are 
similar to those in Fig. 1. The SHG response is found to be enhanced more strongly in the 
presence of both a relaxation-rate gradation and a plasma- frequency gradation (see Fig. 1) 
than for plasmon-frequency gradation alone, especially at low frequencies. As increases, 
the structures in the linear and SHG response both show a shift to lower frequencies. In 
Figs. 1-2, the quantities which can be both positive and negative are plotted in logarithm 
of modulus. When the quantities pass through zero, the logarithm will be very large, thus 
yielding spikes. In addition, we have used the normalized numbers in order to describe a 
general origin of the SHG in metal films rather than a specified metal film. 

The point for achieving the present results is that one needs a sufficiently large gradient 
rather than a crucially particular form of the dielectric function or gradation profiles. Thus, 
it is expected that an enhancement in SHG responses will also be found in compositionally 
graded metal-dielectric composite films in which the fraction of metal component varies 
perpendicular to the film. In the present work, due to the symmetry of the film, we have 
only enhancement for the polarization perpendicular to the film (i.e., parallel to the direction 
of the gradient). In this polarization, the tangential component of electric field E vanishes 
identically. Thus, the continuity of normal component of D [see Eq. (2)] gives rise to 
the enhanced SHG. However, for the polarization parallel to the film, i.e., the tangential 
component does not vanish. In this polarization, it is the continuity of the tangential 
component of E that leads to no enhancement at all [1] . 
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FIGURES 

FIG. 1. (a) Re[e(a;)], (b) Im[e(a;)] (linear optical absorption), (c) Re[x2a)/xi]) (d) Im[x2a;/Xi], 
and (e) Modulus of X2aj/xi versus the normalized incident angular frequency iL)/ujp{d) for the dielec- 
tric function gradation profile [Eq. (1)] with various plasma-frequency gradation profile [Eq. (4)] 
and relaxation-rate gradation profile [Eq. (5)]. Here | ■ • • | denotes the absolute value or modulus 
of • • •. Parameters: 7(00) = 0. 02^^(0) and C-y = 0.0. 



FIG. 2. Same as Fig.l. Parameters: 7(0x3) = 0.02u;p(0) and C^, = 0.6. 
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Fig. 1/Huang and Yu 
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Fig. 2 /Huang and Yu 
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